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Abstract. We employ high-resolution total neutron scattering in conjunction with 
reverse Monte Carlo simulations to examine, in a detailed and unbiased manner, the 
crystal structure of the vacancy-ordered oxide pyrochlore Pb2Ru2060Q 5 in light of 
its structural analogy with proton-ordering in the structures of ice. We find that the 
vacancy and the O' ion are completely ordered, and that the average structure in the 
F43m space group describes the vacancy ordering precisely. We complement these 
results with an examination of the Pb^^ lone pair network using density functional 
electronic structure calculations, and a comparison of the low-temperature lattice-only 
heat capacity of Pb2Ru2060Q 5 with that of other related pyrochlores. 

PACS numbers: 61.05.fm, 61.43.Bn, 
1. Introduction 

In recent years, pyrochlorepj A2B2O6O' compounds (also described compactly as 
A2B2O7) have been explored - in cases where A is a rare-earth ion and B is a transition 
element - for a number of reasons ranging from their unusual magnetic "spin-ice" 
ground states [2j and their anomalous Hall effect behavior [3J to possible realization of 
magnetic monopoles,[l] and because certain pyrochlore compositions, notably Y2lr207, 
can display electron correlation in the presence of spin-orbit coupling. [5] 

When pyrochlores are prepared with lone-pair active ions such as Pb^+ or Bi'^+ on 
their A sites, rather than rare-earth or alkaline earth ions, one can expect interesting 
structural effects. An example of this is the grotesquely complicated low-temperature 
structure of Bi2Sn207 and its numerous polymorphs that form upon heating. P [7] 
Hector and Wiggin[8j reported a low-temperature preparation of Bi2Ti2060' in 2004; 
a compound which prefers to disorder extensively and locally, and remains cubic 
down to a temperature of 2 K rather than undergo a phase transition to a coherently 
disordered ground state. This observation led to the suggestion that just as the 
pyrochlore lattice hinders coherent magnetic ordering [9j so too possibly are coherent 
lone-pair induced displacements hindered in compounds such as Bi2Ti2060' and related 
pyrochlores ITTl [T2] a phenomenon dubbed "charge-ice" [13] in analogy with problems 
of proton ordering in the structure of ice- Jc and ice-//i.|TH |15] The low-temperature 
thermal properties [IS] and detailed scattering studies [T7| have supported the validity 
of this suggestion. 

In contrast to A2B2O6O' pyrochlores, which crystallize in the Fd3m space group and 
whose O' sites describe a diamond lattice, in Pb2Ru2060Q 5 (average structure reported 
from Rietveld refinements of neutron [18] and synchrotron X-ray diffraction [19]) there 
is ordering of O' and a vacancy (depicted by □) so that the compound crystallizes 
in the F43m subgroup of Fd3m, with O' and □ forming a zinc blende lattice. The 
complete ordered structure in space group FASm can be written A2B20i20'n. In 
contrast to A2B2O6O' where the diamond O' lattice has an A ion between every two O', 
in A2B20i20'n, each A ion has an O' and a □ as neighbors, at least in the perfectly 
ordered structure (shghtly further are six O neighbors in the -B2O6 sublattice). 




Figure 1. (Color online) (a) The unit cell of Pb2Ru2060o 5 consisting of corner- 
sharing networks of RuOg octahedra and isolated O'Pb tetrahedra. The ordered 
vacancies (□) are in the 4a(0, 0, 0) positions and are indicated using large white spheres, 
(b) The lone-pairs on the Pb^^ visualized using electron localization function at an 
isosurface of 0.80. For clarity, Ru and O are not displayed in (b). 



The analogy with ice arises in the following manner: the topology of A2O' in the 
pyrochlore structure is analogous to the topology of H2O in ice Jc, with the protons 
in the latter obeying the "2-in, 2-out" Bernal- Fowler rules [H] in a completely random 
way that results in residual entropy. [13, 120] If the title compound were fully ordered, 
the A2(0'n)o.5 network would also be completely ordered with the A ions (Pb^"*") 
choosing to all bond to O' or to □ in a "4-in 0-out" fashion. Hence our use of the 
term "ordered-ice" . Figure[T|a) displays the crystal structure of completely vacancy- 
ordered Pb2Ru2060g 5. The vacancy sites are displayed as large white spheres in order 
to emphasize the Pb2^(0'n)o.5 network in the structure. 

In this contribution, we employ high-resolution neutron scattering tools based on 
analysis using the reverse Monte Carlo method |22] and least-squares refinements 
to establish, in a structurally unbiased manner, the extent of ordering between O' 
and □. We find that the average structural description in the F43m space group of 
completely ordered vacancies is consistent with the results of the reverse Monte Carlo 
simulation of total scattering. An interesting consequence of the ordered vacancies is 
that the stereochemically active lone pairs on the Pb^"*" ions have a location to position 
themselves in a completely coherent fashion as we shall show. We use this study to 
emphasize the importance of employing both the Bragg scattering intensity as well as the 
pair distribution function, G{r) in the reverse Monte Carlo analysis. Finally, we examine 
the lattice part of the low temperature heat capacity for local Einstein modes that have 
been observed in other pyrochlore compounds with incoherent lone pairs. pT| [25] 




Figure 2. (Color online) The neutron TOF Rietveld refinement for Pb2Ru206.5 at 
T = 15K confirms the F43m space group. The data is presented in two panels to 
emphasize the Q-range and the quality of the fit. 



2. Experimental and computational Methods 

The sample for neutron scattering was obtained by careful crushing and grinding of 
Pb2Ru2060Q 5 single crystals, whose average structure by synchrotron X-ray Rietveld 
refinement, and whose electrical transport properties have been previously reported. [TU] 
Time-of-fiight (TOF) total neutron scattering was carried out at the NPDF instrument 
at Los Alamos National Laboratory at 300 K and 15 K. Rietveld refinement was 
performed with the EXPGUI frontend for GSAS.[2l]. Extraction of the PDF G{r) 
with PDFGetN[25] used Qmax = 35 A~^. Reverse Monte Carlo simulations were 
run using RMCProfile version 6 [22] and a 6 x 6 x 6 unit cell with 18144 atoms. 
These simulations were simultaneously constrained by the Bragg profile and the G{r). 
Multiple runs were performed and averaged when necessary to obtain the most accurate 
representation of the fit to data. First principles electronic structure visualization made 
use of the linear muffin-tin orbital method within the atomic sphere approximation 
using version 47C of the Stuttgart TB-LMTO-ASA program. [26] More details of such 
calculations in Pb and Bi-based pyrochlore oxides can be found in reference [13] . 

3. Results and discussion 

The result of Rietveld refinement using neutron TOF scattering is shown in Fig.[2] No 
impurities are present, and the FASm unit cell provides an excellent fit, with details 
into the high-Q range reproduced past Q = 12 A^^. This model contains Pb atoms that 
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Table 1. Average structural parameters for Pb2Ru2060o 5 obtained from Rietveld 
refinement of TOF neutron diffraction data in the space group F43m (No. 216). 
a = 10.25017(2)1/ at 300 K and a = 10.22775(2)1 at 15 K. All sites refine to full 
occupancy. 
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are shifted off the ideal pyrochlore A site into a 16e position at {x^x^x). The O' anions 
occupy Ad sites. Details of the structure obtained from Rietveld refinement of the 300 K 
and 15 K data are displayed in Table[TJ All atoms refine to their full occupancy and all 
the isotropic thermal parameters from the refinement refine to acceptable values. 

With the ideal, average structure, density functional calculations can be performed, 
and used to locate the lone pairs on Pb^+[2Z] using the electron localization function 
(ELF).[2H1 EH] Figure|l|(b) displays the lone pair network, visualized for an isosurface 
value of ELF = 0.80. The lone pairs on Pb^+ are seen to coherently dispose themselves 
in the direction of the □. It is interesting to note that while lone pairs are known to 
possess a volume that is comparable to that of oxide anion, the lone pair to nucleus 
(of Pb^"^) distance is significantly smaller than the distance from Pb^+ to O'. Hence, 
the lone pairs are not actually located in the vacant □ site at 4a(0,0,0). This is in 
keeping with the empirical findings of Galy, Andersson and others. [30] It is instructive 
to compare this ordered lone pair arrangement with what is seen in "charge-ice" 
Bi2Ti2060'. In Bi2Ti2060', which has no vacant site, the lone pairs are obliged to 
displace perpendicular to the O'-A-O' axis, and they do so incoherently, with profound 
structural and thermodynamic consequences. [131 [El E] 

Manually altering the Pb positions or O' ordering in the model reveals that both of 
these values are strongly constrained by the data. In Fig.[3]the effect on the goodness- 
of-fit is plotted as a function of Pb displacement. For an ideal pyrochlore cell, the 
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Figure 3. (Color online) In (a), moving Pb off of its refined position at x — 0.8766 
results in a poorer fit, seen from the poorer fit (larger chi^) for data acquired at 15 K 
and 300 K. In (b), O' occupancy perfectly ordered (left) and increased to simulate 
disorder spilling onto the vacant sites. Again, even small amounts of disordered O' 
give a poorer Rietveld refinement. 



Pb would lie at the left-most end of the plot with x = | = 0.875. Refinements at T = 
15 K and 300 K both suggest Pb positions near 0.877, which are about 0.04 and 0.03 A 
displaced, respectively, off their ideal sites and away from the nearby O'. 

Moving fractional O' occupancy onto the 4a vacant □ site in our model is equivalent 
to introducing site disorder, where having full occupancy on both sites would mimic the 
O' sublattice of the full pyrochlore structure. The effect of increasing O' occupancy 
on the □ site (with corresponding decrease at the O' Ad position) is seen in Fig.[3]^b). 
Again, any occupancy on 4a results in an increased x^, implying that Bragg scattering 
supports O' being fully ordered, even at 300 K. This strong proclivity to order the O' 
and, in turn, the Pb displacements, is seen in related solid solutions Ph2-xLnxIin20i^y 
{Ln = Nd, Gd) which can tolerate only a small amount {x < 0.2) of A-site substitution 
before the disorder drives formation of a separate Fd3m phase. [31] 

The quality of the fit to the Rietveld-refined structure and the stabilities of the Pb 
position and O' ordering provide evidence that Pb2Ru2060o 5 appears structurally to 
be consistent with the unit cell, unlike the full-0' pyrochlores such as Bi2Ti2060' where 
the difficulty in accomodating the lone pairs results in incoherent structural disorder 
and off-centering. In Bi2Ti2060', local structure analysis reveals details of the short- 
range structure that are difficult to determine from average structure probes that rely on 
Bragg scattering alone. In a system such as Pb2Ru2060o 5 where structural frustration 
should be removed by O' ordering, it is of interest to ask whether the local structure 
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Figure 4. (Color online) Reverse Monte Carlo (RMC) fits to the Bragg profile (a) 
and G(r) (b) obtained from neutron total scattering at 15 K. Pair correlations are 
extracted from these RMC supercells. 

exhibits any hallmarks of a disordering that is not seen in the Bragg scattering/Rietveld 
analysis. 

Reverse Monte Carlo fits to the neutron TOF total scattering are shown in Fig.|4j 
The Bragg profile retains the long-range periodicity and average atomic displacement 
parameters of the structure, while the G{r) constrains pairwise distances and reproduces 
the local structure. By simultaneously fitting both datasets, we are left with a large 
supercell that provides a snapshot of the material. Structural details of this model 
can be analyzed statistically to determine the precise local tendencies supported by the 
experimental data. 

The partial pair distribution functions in Fig.|5](a) display the Pb-0' distances in 
the unit cell npb^o'i^) the distances of Pb to the vacant site npb-n{f^)- The shift 
versus r shows that Pb-0' distances are longer, so Pb nuclei have displaced off their 
ideal site toward the vacant □ position. 

Partial pair distribution functions can also reveal ordering on the O' sublattice 
when vacancy "atoms" with zero scattering cross-section are placed on the vacant O' 
positions. The vacancies are freely allowed to swap with O' during the RMC simulation. 
For a simultaneous fit to Bragg peaks and the G{r), ordering is retained: the r = 4.45 A 
nearest-neighbor distance between occupied and vacant O' sites has no intensity for the 
no'-o'{i^) distribution in Fig.jsj^b), indicating that no two O' atoms share neighboring 
sites. However, if the RMC simulation is performed as a fit to the G{r) only (and the 
Bragg profile is ignored) then the O' ordering begins to melt and O' can be found on the 




Figure 5. (Color online) In (a), pair correlations from the RMC supercell show the 
Pb-0' and Pb-vacancy distances. The peaks are offset, indicating a shift of Pb toward 
O'. In (b), O'-O' and O'-vacancy pair correlations show that O' are ordered when 
RMC simulations are constrained by the G(r) and Bragg profile. A fit unconstrained 
by Bragg scattering (c) becomes partially disordered on the O'-vacancy sublattice-the 
G(r) is insufhcient to reproduce the ordered distribution via RMC. 



vacancy sites, as seen in Fig[5[c). Both models are, at some level, "correct" because they 
reproduce the respective experimental data. However, the data guiding the G(r)-only fit 
is incomplete. The O' ordering prescribed by the Bragg pattern and verified in Fig.|3](b) 
is unambiguous, so the fully ordered case is confirmed. Because RMC is a stochastic 
technique, it reveals that the G{r) alone does not contain sufficiently robust information 
about the O' occupancy to prevent O' from spilling onto disordered sites. Only when 
the additional Bragg constraint is applied is the fully ordered model reproduced from 
the data. 

Finally, it is of interest to see whether the ordered arrangement of Pb^"*" ions and 
their associated lone pairs result in distinctly different low-temperature heat capacity 
behavior in Pb2Ru2060o 5 when compared to a compound with disordered lone pairs 
such as Bi2Ti2060'. We display such a comparison in Fig.[6] The scheme used to 
display the scaled heat capacity, Cp/T^ vs. T allows local Einstein-like modes to be 
disinguished. [HH] and the temperature of the peak(s) can be read off as energies of the 
Einstein mode. Liu and Lohneysen|M] have suggested a general scaling law wherein the 
lower the temperature at which Cp/T^ displays a maximum, the larger the peak height 
as well. They have made such correlations for a number of amorphous and crystalline 
systems and Bi2Ti2060' (with its hightly disordered Bi20' network p!7]) displays, of 
all studied crystalline systems, one of the lowest peak temperatures and largest peak 
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Figure 6. (Color online) Lattice-only heat capacity, with the electronic contribution 
subtracted, of Pb2Ru2060Q5, displayed in a manner that emphasizes Einstein-like 
features at low temperatures. For comparison, the similarly scaled heat capacity of 
"charge-ce" Bi2Ti2060' (from reference [TB|) and of pyrochlore Y2Ti207 (taken from 
reference [32]) which does not have lone pairs, is also displayed. 

amplitudes. In contrast, we see from Fig. [6] that in Pb2Ru2060Q 5, (using published 
data from reference [12] after subtracting the electronic 7T contribution) the peak 
occurs at significantly higher temperatures and is significantly smaller in amplitude. 
It is nonetheless larger than what is seen for Y2Ti207[32] which has no lone pairs. 
The question is why, in a completely ordered structure, is there an Einstein-like mode 
at all? The suggestion made recently by Safarik et a/. [35] is that Einstein-like modes 
appear even in crystalline materials when van Hove peaks in the phonon densities of 
state cross the chemical potential. In any event, the trend seen in Fig.|6] is consistent 
with the finding that Pb2Ru2060Q5 has a well-ordered network of A-site atoms unlike 

Bi2Ti2060'. 

4. Conclusions 

Our study of the ordered-ice pyrochlore Pb2Ru2060o5 has found that Rietveld 
refinements dictate O'-D occupancy to be fully ordered, and Pb^"*" is displaced toward 
O' (0.04 A at T = 15 K). This ordering implies that Pb2Ru2060Q 5 does not possess 
the same geometric frustration seen in fully-occupied A2B2O7 pyrochlores with lone- 
pair active cations on the A sites. To compare Pb2Ru2060o5 to those systems, 
where local structure analysis provide a view of incoherent distortions, we performed 
RMC simulations which revealed the presence of coherent Pb^-|- off-centering only. 
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Simulataneous fits to the Bragg intensity and G{r) are required to reproduce the 
vacancy ordering, which is not strongly constrained by the G{r) alone. First-principles 
calculations and heat capacity measurements corroborate to show that the O'-D ordering 
provides an outlet for (and perhaps is constrained by) the pointing of lone pairs into 
vacant 4a positions. 
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